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Abstract: A large-scale afforestation project has been carried out since 1999 in the Loess Plateau of 
China. However, vegetation-induced changes in land surface temperature (LST) through the changing land 
surface energy balance have not been well documented. Using satellite measurements, this study quantified 
the contribution of vegetation restoration to the changes in summer LST and analyzed the effects of 
different vegetation restoration patterns on LST during both daytime and nighttime. The results show that 
the average daytime LST decreased by 4.3°C in the vegetation restoration area while the average nighttime 
LST increased by 1.4°C. The contributions of the vegetation restoration project to the changes in daytime 
LST and nighttime LST are 58% and 60%, respectively, which are far greater than the impact of climate 
change. The vegetation restoration pattern of cropland (CR) converting into artificial forest (AF) has a 
cooling effect during daytime and a warming effect at nighttime, while the conversion of CR to grassland 
has an opposite effect compared with the conversion of CR to AF. Our results indicate that increasing 
evapotranspiration caused by the vegetation restoration on the Loess Plateau is the controlling factor of 
daytime LST change, while the nighttime LST change is affected by soil humidity and air humidity. 
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1 Introduction 


Arid and semi-arid regions cover 41% of the Earth's land surface, with 38% of the human 
population living in the regions. These regions are more ecologically vulnerable and sensitive to 
climate change and human activity, which can affect the health of humans (Reynolds et al., 2007; 
Yang et al., 2015, 2016). A series of large-scale ecological restoration projects have been carried 
out throughout the world to address potentially serious environmental and ecological issues, and 
China has also taken efforts in recent decades on the issue (Yi et al., 2009; Allan et al., 2013). In 
1999, a large-scale vegetation restoration project, Grain for Green (GFG) which is the largest 
eco-environmental construction project and most important nature conservation program in arid 
and semi-arid regions of China (Liu et al., 2008; Feng et al., 2013), was launched by the Chinese 
government. GFG involves an effort to return croplands on steep slopes to forest lands and 
grasslands by using a government payment scheme, which engages millions of rural households 
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as core agents of project implementation (Lii et al., 2012; Feng et al., 2013). As the core region of 
GFG, the Loess Plateau is an ecologically vulnerable area, sensitive to climate factors (e.g., 
periodic heavy rainfall) and intensified human activities (e.g., over-grazing and coal mining) 
(Zhang et al., 2008; Liang et al., 2013). However, the vegetation coverage of the Loess Plateau 
has gradually improved since the implementation of GFG (Xin et al., 2008; Zhou et al., 2009), 
and over 6x10° km? of croplands have been converted into forest lands and grasslands during 
1999-2010 (Li et al., 2016). 

Numerous studies have attempted to evaluate the impact of the vegetation restoration project 
on the Loess Plateau from different perspectives. For example, Gao et al. (2017) used the Budyko 
framework to simulate the water cycle on the Loess Plateau during 1990-2014. With an increase 
in vegetation coverage, both regional precipitation and actual evapotranspiration showed a 
significant increasing trend, and the increasing rate of precipitation was greater than 
evapotranspiration. Deng et al. (2016) found that the storage of soil carbon increased by about 
21.3% in the topsoil (0-20 cm) of the re-vegetated lands. Zhai et al. (2015) investigated the 
relationship between vegetation restoration and land surface albedo on the Loess Plateau for the 
period 2000-2010 and believe that human-induced land use change is the control factor for 
albedo reduction. In addition, some studies discussed the effects of vegetation change on energy 
and material recycle from different perspectives, including energy budgets (Chapin et al., 2005; 
Mevicar et al., 2007), soil organic matter (Xiao et al., 2017), soil humidity (Jian et al., 2015), and 
soil erosion (Zhou et al., 2006; Zhang et al., 2016). Based on these studies, we conclude that the 
vegetation restoration project has greatly affected the environment and material recycling 
processes on the Loess Plateau. However, the changes in land surface temperature (a dominant 
factor in fragile ecological regions) caused by the increased vegetation coverage has not been 
well documented. 

As the direct driving force in the exchange of long-wave radiation and turbulent heat fluxes at 
the surface-atmosphere interface, the land surface temperature (LST) is one of the most important 
parameters in the physical processes of surface energy and water balance at scales from local to 
global (Anderson et al., 2008; Karnieli et al., 2010). On the one hand, the changes in LST will 
directly affect soil moisture content and soil evaporation, thereby indirectly affecting regional 
vegetation growth, agricultural production, and the water cycle (Zhang et al., 2014; Harris et al., 
2017). On the other hand, the LST is an important indicator for predicting drought, climate 
change, and earthquakes (Choudhury et al., 2007; Amiri et al., 2009; Vancutsem et al., 2010). The 
LST is the joint outcome of sensible heat and latent heat, which is influenced by albedo, 
evapotranspiration, surface roughness, and other factors (Xue et al., 2017). According to Peng et 
al. (2014), vegetation restoration changes the albedo and thus changes the amount of absorbed 
solar radiation at the ground surface. At the same time, vegetation restoration generally increases 
evapotranspiration (ET) and causes cooling effect (Jackson et al., 2008). Therefore, vegetation 
restoration may result in both cooling effect and warming effect, which is related to the vegetation 
restoration pattern and latitude (Lee et al., 2011). The biophysical effects of vegetation restoration 
on local climate can be much larger than the small global cooling effect resulting from the uptake 
of CO2 by growing vegetation (Lee et al., 2011). However, these biophysical effects are also 
complex and depend on "background" climate (Wang et al., 2011). In general, although the 
large-scale vegetation restoration project has been conducted for more than 15 years, there is a 
complex relationship between vegetation restoration and LST, and as such, the effects of 
vegetation change on the LST on the Loess Plateau are not clear and need to be further elucidated. 

In this study, we investigated how vegetation restoration affects the LST on the Loess Plateau, 
using satellite-derived LST datasets from MODIS instruments in summer (June to August) during 
2000 to 2013. These LST data depend on the radiative properties of the land surface (Wan, 2008; 
Zhou et al., 2012) and, therefore, have a larger diurnal amplitude than the observation data from 
meteorological stations (Jin and Dickinson, 2010). The objects of this study are: (1) to reveal the 
characteristics and controlling factors of the LST change in the vegetation restoration area; and 
(2) to quantify the effects of different vegetation restoration patterns on the LST during both 
daytime and nighttime. 
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2 Materials and methods 


2.1 Study area 


The Jinghe River Basin, a small region with an area of 4.54x10* km’, was chosen as the study 
area to minimize the effect of altitude and precipitation on the LST. The basin is located in the 
southern Loess Plateau of China (34°14'—-38°10'N, 105°49’-108°58’E; Fig. 1), dominated by a 
semi-arid continental monsoon climate. During the period 1980-2015, the annual summer 
precipitation varied from 190 mm in the northwest to 250 mm in the southeast part of the basin, 
with a mean value of 225 (+30) mm. It has a typical water-limited landscape, with actual 
evapotranspiration accounting for 85% of the precipitation; this percentage increases along the 
precipitation gradient (Jin and Dickinson, 2010). The main soil types in the basin are Lossiah soil 
and dark Lossiah soil, which are thick soil layers (50—80 m), loose in structure, easily collapsed 
and eroded. The land-use types are primarily cropland, grassland, forest, and sandy land. In 
Jinghe River Basin, a large number of ecological restoration projects has been implemented in the 
past 30 years since the 1980s. However, these costly efforts yielded little success during the first 
20 years (Peng et al., 2015). In the past 10 years, vegetation restoration has begun to achieve 
certain effects, and the average summer NDVI has increased from 0.4 in 2000 to 0.7 in 2013. 
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Fig. 1 Location map and change in NDVI in the study area. (a) NDVI in 2000, (b) NDVI in 2013, and (c) NDVI 
change from 2000 to 2013. 


2.2 Data 


We used the MODIS collection-3 products downloaded from NASA (https://www.nasa.gov/topics/ 
earth/index.html) for the LST, surface albedo, and ET. These products have been used extensively 
in a variety of areas and are proven to be of high quality (Wan, 2008; Zhou et al., 2012). For LST, 
we used the 8-day average LST from Terra-MODIS (MOD11A2) products in summer (June to 
August) from 2000 to 2013. The LST data cover both daytime (10:30 am) and nighttime (22:30 
pm). In this study, we corrected the noise resulting from cloud contamination and topographic 
differences, and only pixels with LST error <1 K were used. The MODIS surface albedo products 
(MCD43A3) include white sky albedo and black sky albedo, with a 0.5-km spatial resolution and a 
16-day time resolution from 2000 to 2013. We used the white sky albedo shortwave broadband in 
this study, and the bias of the MODIS albedo was mostly less than 5% (Liu et al., 2009). The 
MODIS ET data (MOD16A2) uses land cover, albedo, air temperature, radiation, humidity, and 
leaf area index as input data, with a 0.5-km spatial resolution and an 8-day interval. The mean bias 
of the MODIS ET is 0.3 mm/d compared to ET observations from the eddy flux tower (Mu et al., 
2011). 

The NDVI data and land cover map (2000, 2005, 2010 and 2013) at 1-km resolution were 
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downloaded from Resource and Environment Data Platform  (http://www.resdc.cn/ 
DataList.aspx), of which NDVI data come from vegetation sensor with monthly interval and land 
cover data come from Landsat TM and Landsat ETM. In this study, we adopted the China land 
resource classification system for the land cover data (Zhang et al., 2010). 

Meteorological data, such as relative air humidity, observed land surface temperature, and other 
data needed to calculate the net radiation were collected from six national weather stations. The 
relative soil humidity was collected from 57 national agricultural meteorological stations in the 
study area from 2000 to 2013. Meteorological data and relative soil humidity data were 
downloaded from Meteorological Data Center of China Meteorological Administration 
(http://data.cma.cn/). 


2.3 Methods 


In the study area, we defined the region where the average NDVI increased by more than 0.2 
during 2000 to 2013 as the vegetation restoration area (VRA) (blue regions in Fig. 1c) (Cao et al., 
2017) and the area where the NDVI change was less than 0.2 as the non-afforestation area (NAA) 
(yellow regions in Fig. 1c). In addition, irrigated areas, mountains (slope>15°; Deng et al., 2016), 
and urban areas were excluded (the areas of VRA and NAA account for 52% and 46% of the 
basin area, respectively). The average elevation difference between VRA and NAA is 178 m, less 
than the calculation limit of 200 m (Peng et al., 2014). 

Because the LST in VRA was under the influence of vegetation restoration and climate, the 
LST change in NAA without the influence of vegetation was set as LST climate, the LST in VRA as 
LSTvyra, and the effect of vegetation on LST as LST vegetation. The relationship between LSTvra 
and LST vegetation is Shown as Equation 1. 


LS Tyra = ESE aP ESE Es (1) 


Remote LST was defined as the average of the daytime LST and nighttime LST. All daily 
meteorological data were calculated to obtain a monthly average, which was finally calculated to 
obtain the seasonal average data for each year. The changing trend in the time series data was 
analyzed with the Mann-Kendall method. 


3 Results 


3.1 LST change in VRA 


The LST in VRA significantly decreased during the summer and exhibited a clear correlation with 
NDVI (correlation coefficient=0.7). We can clearly see in Figure 2a that the average NDVI in 
VRA increased from 0.40 in 2003 to 0.67 in 2013. With the increase in NDVI, summer LST in 
VRA decreased from 25.2°C in 2000 to 23.5°C in 2013. LST is the average value of daytime LST 
and nighttime LST, and these two factors together determine the changes in LST. The 
Mann-Kendall trend analysis indicates that the daytime LST significantly decreased from 2000 to 
2013 (Test Z= —2.29, 95% confidence level), and the average temperature decrease was 0.3°C/a. 
The LST nighttime slightly increased during the same period (Test Z=1.20, did not pass the 
significance test of 95%), with an average rise of 0.1°C/a. The daytime cooling effect is 
significantly higher than the nighttime warming effect in VRA, meaning that the significant 
decreasing trend in LST from 2000 to 2013 is mainly attributed to the daytime LST change. We 
can also see from Figures 2b and c that the daytime LST and nighttime LST in VRA both showed 
large fluctuations from 2000 to 2013, which may be attributed to climate change and the low time 
resolution (8-day) of the LST product. 


3.2 Effects of vegetation restoration on LST 


The LST change in VRA is affected by vegetation restoration and climate factors. To reveal the 
impact of the vegetation change on the LST, we assume that the LST in NAA is only affected by 
climatic factors, and LST imate is used to reflect the influence of climatic factors on the LST of 
VRA. Based on remote LST data, we plotted the inter-annual variation curve of LSTvyra, 
LST climate, and LST vegetation in the summer from 2000 to 2013 (Figs. 3a and b). 
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Fig. 2 Summer land surface temperature (LST) change during 2000 to 2013. (a) LST and NDVI change, (b) 
daytime LST change, and (c) nighttime LST change. 


In daytime, LSTclimate shows a slight downward trend (Test Z= —0.99, did not pass the 
significance test of 95%), with an average decrease of 0.124°C/a. However, LSTvra decreased 
significantly during the same period, with an average decrease of 0.311°C/a, which is far greater 
than the annual decrease trend of LST climate. These results suggest that the significant decreasing 
trend in daytime LSTvra could be mainly attributed to vegetation restoration, and this cooling 
effect became stronger with the increase in NDVI (Fig. 3a). Based on Equation 1, LST vegetation 
experienced an extremely significant increase from —3.99°C in 2000 to —6.13°C in 2013 (Test 
Z=3.28, 99% confidence level), and the cumulative cooling effect exceeded 2.5°C in summer 
from 2000 to 2013 (Fig. 3a). By comparing the effect of LST vegetation and LST climate on LSTvra, we 
conclude that vegetation restoration contributed 58% of the decrease in daytime LSTvra, which 
was far greater than the contribution by the climate change (42%). 

In nighttime, both LSTvra and LST climate show a slight upward trend from 2000 to 2013, and 
the upward trend in LSTyra is greater than that in LSTcimate. This result means that the change in 
vegetation helps increase the nighttime LST in VRA. Based on Equation 1, vegetation change 
contributes a little to the warming effect at nighttime, with an average increase of 0.06°C/a during 
2000-2013 (Fig. 3b). Based on the trend in LST vegetation, LST climate, and LSTvra, we also can 
conclude that vegetation restoration contributed 60% of the increase in nighttime LSTvyra. 
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Fig. 3 Summer LSTvra, LSTvegetation, and LSTciimate change in vegetation restoration area during 2000-2013. (a) 
daytime and (b) nighttime. 
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3.3 Effects of different land cover types on LST 


Converting slope cropland (CR) into artificial forest (AF) and grassland (GR) are the main 
patterns of vegetation restoration on the Loess Plateau. This study used remote LST data and land 
cover data to analyze the impact of CR, AF, and GR on daytime and nighttime LST in VRA (Fig. 
4). The average LST of AF was significantly lower than those of CR and GR during the daytime, 
at about 2.2°C between AF and CR and 2.4°C between AF and GR. The difference in daytime 
LST between CR and GR is small; however, the LST of GR is always higher than that of CR, 
with an average gap of 0.21°C. During nighttime, although the differences in LST among AF, CR, 
and GR are only 0.07°C—0.63°C, the nighttime LST of AF is always slightly higher than those of 
CR and GR, with differences of 0.56°C and 0.63°C, respectively. The land surface temperature of 
different land cover is affected by both climate and vegetation factors. According to Equation 1, 
the average LST of AF is —5.6°C, which is significantly smaller than those of GR (—3.2°C) and 
CR (-3.4°C). This indicates that AF has a stronger cooling effect than CR and GR in daytime. 
Although AF, CR, and GR all show a small warming effect at nighttime, the warming effect of AF 
(0.70°C) is always greater than those of CR (0.12°C) and GR (0.05°C). Compared with the other 
years, the LST is significantly varied under different land covers during daytime and nighttime in 
2005. This is because the LST in 2005 had larger fluctuations than in the other years, which is 
mainly attributed to climate factors (Fig. 3). 

The LST changes of daytime and nighttime in different vegetation restoration patterns were 
analyzed in this study (Fig. 5). For CR-AF (conversion from CR to AF), the cooling effect was 
observed in the daytime and the warming effect occurred at nighttime, which is consistent with 
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Fig. 4 Effects of croplands, artificial forest, and grasslands on LST from 2000 to 2013. (a) daytime and (b) 
nighttime. CR, cropland; GR, grassland; AF, artificial forest. 
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Fig. 5 LST changes of daytime and nighttime in different vegetation restoration patterns 


the result of the change in LST in VRA. However, for CR-GR (conversion from CR to GR), the 
diurnal change in LST has the opposite effect compared with CR-AF. In daytime, the cooling 
effect of LST of CR-AF is 2.2°C. In contrast, the average annual nighttime LST between CR and 
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AF was 0.2°C, indicating that the vegetation restoration of CR-AF warms the land surface during 
nighttime but with a smaller magnitude than the daytime cooling effect. The warming and cooling 
effect of GR are always lower than those of CR (Fig. 4), so the vegetation restoration of CR-GR 
increase LST during the day (0.7°C) and decrease LST at nighttime (0.1°C). 


4 Discussion 


4.1 Influence factors of LST change in daytime and nighttime 


Vegetation restoration has the cooling effect in daytime and the warming effect at nighttime in 
VRA (Fig. 2). The asymmetric diurnal cycle of LST results from different energy balance 
processes (Dai et al., 1999; Zhou et al., 2007). LST is determined by sensible heat (Allan et al., 
1998); therefore, the LST can be expressed as a function of sensible heat (H) in the study area 
(Eq. 2). 


LST = f (AR). (2) 
Based on the surface energy balance equation, the H is given as Equation 3. 
H=R,-LxET-G, (3) 


where H is the sensible heat flux (MJ/(m?-d)); Rn is the net radiation (MJ/(m2-d)); ET (mm) and G 
(MJ/(m?.d)) are the evapotranspiration and soil heat flux, respectively, and L is the latent heat of 
vaporization of water, which is a constant (L is about 2.45 MJ/kg) (Allan et al., 1998). 
Rn was calculated using the Penman-Monteith equation: 
R, =Rs -Ru (4) 
where Ras and Ru are net short-wave radiation (MJ/(m?-d)) and net long-wave radiation 
(MJ/(m?.-d)), respectively. Rns and Rn can be given as Equations 5 and 6 (Allan et al., 1998). 


Rys z (l-a@)(a, aD “R, ? (5) 
N 
bn 
1.35(a, +) 
oT. +0T.. s N 
R =—™ = (0.34-0.14 x Je, )(————— 0.35) . 6 
nl 2 ( mS e, )( +b, ) ( ) 


The meanings and values of each coefficient in Equations 5 and 6 can be found in Liu and 
McVicar (2012). The albedo (a) can be set as a different value for different NDVI (Zhai et al., 
2015), and changes in the albedo will further affect R, and LST (Peng et al., 2014). Therefore, 
based on MODIS albedo data, we used the annual summer surface albedo from 2000 to 2013 to 
reflect the effect of vegetation change on LST. 

The daytime LST is mainly controlled by R, and LxET (Eq. 3; R, and LxET are much larger 
than G in daytime and this study ignores the effect of G on sensible heat flux) (Allan et al., 1998). 
Summer extra-terrestrial solar radiation is assumed to be similar in VRA during 2000-2013 
(17.94-21.91 MJ/(m?-d)). Hence, the surface albedo determines the amount of net radiation under 
the same climate conditions (Eqs. 5 and 6). According to Zhai et al. (2015), there is a significant 
negative correlation between NDVI and albedo; thus, vegetation restoration will increase the R, 
and bring warming effect in daytime (Eqs. 5 and 6). During the same time, the increase in 
vegetation will also increase the evapotranspiration and bring a cooling effect (Zhai et al., 2015). 
So, the effect of vegetation change on the daytime LST is determined by both the warming 
(caused by albedo) and cooling effect (caused by ET). Based on MODIS data, this study analyzed 
the change in albedo and ET in the VRA from 2000 to 2013. With the increase in vegetation, the 
ET in summer showed an extremely significant increase (Z=4.38, 99% confidence level), and the 
annual increase trend is 0.1 mm/d (Fig. 6a). During the same period, the albedo significantly 
decreased from 0.193 to 0.152 (Z= —4.05, 99% confidence level), a decrease of 21% and an 
average annual decrease of 0.0023 (Fig. 6b). The effects of albedo and ET change on heat flux are 
calculated based on Equations 3—6 (Fig. 7). We can clearly see in Figure 7 that vegetation 
restoration led to an increase in both net radiation and latent heat in VRA; however, the increase 
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trend in latent heat is significantly greater than the net radiation, and the gap between latent heat 
and net radiation gradually expands, indicating that the cooling effect caused by the ET increase 
is greater than the warming effect caused by the albedo decrease. Thus, ET is the control factor of 
the daytime LST change in VRA. 

In general, land surfaces absorb and store energy from the atmosphere during the day and 
release energy during the nighttime. ET is negligible in nighttime and LST is closely related to the 
energy stored during the daytime and the near-surface atmosphere boundary layer (Wan, 2008; 
Lee et al., 2011). Based on observation data from national weather stations and agricultural 
meteorological stations, the change in relative soil humidity and relative air humidity during 2000 
to 2013 was analyzed. During the period, the relative soil humidity in VRA increased from 58% 
to 67%, with an annual average increase of 0.48% (Fig. 8a). Increased soil humidity resulted in 
more daytime heat storage and more nighttime heating. The increase in air humidity near the 
surface means more downward long-wave radiation received from the atmosphere and a reduction 
in the outgoing long-wave radiation from the land surface (Fig. 8b). This long-wave radiative 
imbalance has a stronger effect during nighttime when the boundary layer is thinner and more 
stable (Zhou et al., 2009). 
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Fig. 6 The change in ET(a), albedo (b), relative soil humidity (c), and relative air humidity (d) in the vegetation 
restoration area during 2000 to 2013. 
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Fig. 7 Effect of ET and albedo change on the heat flux during 2001-2013 (fixed albedo and ET in 2000 as 
initial state, the change in ET and albedo means subtract the initial value from 2001 to 2013, respectively) 
4.2 Influence factors of vegetation restoration patterns on LST 
Different afforestation patterns have different effects on ET and albedo (Liu et al., 2009). The ET 
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Fig. 8 The change in relative soil humidity (a), and relative air humidity (b) in the vegetation restoration area 
during 2000 to 2013. 


of AF is indeed higher than that of CR in 50% of the VRA and higher than that of GR in 68% of 
the VRA, while the albedo is lower in 8% of CR and 9% of GR, respectively (Fig. 9). This clearly 
indicates that AF across the VRA absorbs more incoming radiation and dissipates more energy as 
latent heat. Thus, the artificial forest induces cooling during daytime. The surface incoming solar 
radiation in summer is between 17.94 and 21.91 MJ/(m?°-d) on the Loess Plateau. The average 
difference in albedo is 0.014 between CR and AF, and 0.003 between GR and CR. Hence, the 
extra solar energy absorbed by the change in vegetation cover from CR to AF (CR-AF) is 0.25- 
0.31 MJ/(m-d), while the solar energy dissipated by the change from CR to GR (CR-GR) is —0.05 
to —0.07 MJ/(m?.d). The ET change resulting from these two vegetation restoration patterns is 
0.72 and —0.16 mm/d respectively, which is equal to an extra energy dissipation through latent 
heat of about 1.76 and —0.39 MJ/(m?.d), respectively. So, based on Equations 3—6, vegetation 
restoration patterns of CR-AF will result in a cooling effect during the day due to more energy 
dissipation, while CR-GR will result in a warming effect due to more energy absorption. 
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Fig.9 Average ET and albedo of CR, AF, and GR during 2000-2013 


The impacts of different land use types on soil humidity and air humidity are different in 
nighttime. According to Xu et al. (2008), under the same climatic conditions, the relative soil 
humidity of AF in the study area is about equal to CR, which is 1.5 to 2.0 times that of GR. The 
relative air humidity (150 cm above vegetation) of AF is 13% more than CR in summer, while the 
relative air humidity of CR and GR is similar (Li et al., 2010). This result indicates that CR-AF 
will increase the ability to absorb radiation (increase in air humidity), while CR-GR will decrease 
the ability to store heat (decrease in soil humidity). Thus, the vegetation restoration of CR-AF 
increases LST in nighttime, while CR-GR reduces LST. In addition, the nighttime warming effect 
of CR-AF could be magnified as a result of reduced atmosphere turbulence from a more stable 
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stratification over vegetation (Lee et al., 2011). This reduces the heat dissipation from AF cover 
types compared with the open lands of CR and GR. 


4.3 Uncertainty analysis 


The effect of vegetation on LST between the Jinghe River Basin and other regions is compared to 
evaluate the accuracy of the calculated results (Table 1). Table 1 shows that our results match well 
with Hao et al., (2016), Peng et al. (2014) and Potchter et al. (2008). However, in Tibetan Plateau, 
vegetation greening increases the local LST due to the huge decrease in albedo (Shen et al., 
2015). In general, effects of vegetation cover on LST are varied with the vegetation restoration 
pattern and latitude (Lee et al., 2011). 


Table 1 Comparison between results from other studies 


Location Period Vegetation type Average LST change (°C) Source 
Jinghe River Basin 2000-2013 Restoration —1.65 This study 
Tarim River Basin 2001-2013 Natural -3.02 Hao et al. (2016) 

China 2003-2010 Restoration —0.90 Peng et al. (2014) 
Southern Israel 2004-2006 Natural —2.00 Potchter et al. (2008) 
Tibetan Plateau 2000-2010 Natural +0.78 Shen et al. (2015) 


To test the rationality of LST of MODIS data, the study analyzed the remote LST data and 
observed LST data from six national weather stations in VRA during 2000 to 2013. The remote 
sensing data has a good linear correlation with the observed data (Fig. 10a). Although, the 
MODIS data is usually 0.68°C—1.02°C lower than the actual LST, they both exhibit a similar 
trend. Therefore, the MODIS data in VRA can meet the requirements of the study. 
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Fig. 10 Scatterplot of uncertainty analysis. (a) Remote LST and observed LST and (b) estimated Rn and 
observed Rn. 


In our research, we obtained the summer R, from the Penman-Monteith equation. The observed 
data from the weather station (there is only one weather station in VRA that can observe the net 
radiation and observation time start from 2006) was analyzed to test the rationality of the 
estimated Rn. Figure 10b shows that the correlation coefficient between the estimated R, and the 
observed R, reached 0.81, and the error between observed data and estimated data is 4%-9%. 
Therefore, the estimated R, can meet the requirement of the study. 

Although many factors are considered, this study still has limitations about the effects of 
vegetation changes on LST. For example, vegetation greening can lead to the increase in surface 
roughness and more solar energy could be absorbed by the surface accordingly. Our results may 
overestimate the cooling effect caused by vegetation change. This work also ignores the effect of 
human activities on LST, such as irrigation. In addition, it must be noted that the Loess Plateau is 
a very large and geodiverse region; thus, the research results in Jinghe River Basin may not apply 
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to the entire Loess Plateau. 

Vegetation restoration in Jinghe River Basin has gradually reduced the LST in summer. 
Because the temperature of the deep soil is relatively stable, the decrease in LST will directly 
reduce the soil temperature gap, thus, reducing soil evaporation and increasing the soil moisture 
content (Kidron and Kronenfeld, 2016; Ren et al., 2017). Increased soil moisture caused by 
vegetation change is an important reason for reduced runoff on the Loess Plateau (Fu et al., 
2003); and this could lead to the uncertainty in regional water resources management as well 
(Feng et al., 2016). 


5 Conclusions 


Based on the remote LST data and land cover data, we concluded that the cooling effect and 
warming effect of PF, CR and GR are gradually decreasing. Therefore, the vegetation restoration 
pattern of CR-PF resulted in the cooling effect in the daytime and the warming effect in 
nighttime, while the vegetation restoration pattern of CR-GR has the opposite effect. 

Based on surface energy balance equation, this paper discussed the effect of vegetation 
restoration on energy budget. With vegetation change, net radiation (inflected by albedo) and 
latent heat (inflected by ET) in VRA show an increasing trend; however, the average increase 
trend of latent heat (—2.4 MJ/(m?-d) was significantly greater than the net radiation (0.09 
MJ/(m?.d)). Thus, the daytime LST in VRA was decreased. For vegetation restoration pattern of 
CR-AF and CR-GR, CR converted into AF means that about 1.5 MJ/(m?-d) of energy is released, 
while CR converted into GR means that about 0.3 MJ/(m?.d) of energy is absorbed. Therefore, 
vegetation restoration pattern of CR-AF decreased LST during the day, while CR-GR increased 
LST. The increase of nighttime LST in VRA can be attributed to the increase of relative soil 
humidity and relative air humidity. 

Although many factors are considered in this study for LST attribution analysis, there are still 
some limitations. It should be noted that the Loess Plateau is a very large and geodiverse region, 
as such, the research results in Jinghe River Basin may not apply to the entire Loess Plateau. 
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